Giant Charge Relaxation Resistance in the Anderson Model 
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We investigate the dynamical charge response of the Anderson model viewed as a quantum RC 
circuit. Applying a low-energy effective Fermi liquid theory, a generalized Korringa-Shiba formula 
is derived at zero temperature, and the charge relaxation resistance is expressed solely in terms of 
static susceptibilities which are accessible by Bethe ansatz. We identify a giant charge relaxation 
resistance at intermediate magnetic fields related to the destruction of the Kondo singlet. The 
scaling properties of this peak are computed analytically in the Kondo regime. We also show that 
the resistance peak fades away at the particle-hole symmetric point. 
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In recent years, an experimental effort has been de- 
voted to manipulate and measure electrons in nanocon- 
ductors in real time [T] . At frequencies in the GHz range 
and cryogenic temperatures, current and noise measure- 
ments provide information on the quantum dynamics of 
charge carriers. Experiments have followed essentially 
two directions, by either using on-chip quantum detec- 
tors [2] or by directly measuring the current using low 
noise amplifiers [3]. In an original experiment, Gabelli 
et al. have realized the quantum equivalent of an 
RC circuit with a quantum dot connected to a spin- 
polarized single-lead reservoir and capacitively coupled 
to a metallic top gate, this setup being later used as 
a single-electron source [5]. By applying an AC mod- 
ulation to the gate voltage, they measured the admit- 
tance of the dot at low frequency. A comparison with 
the classical RC circuit allows to extract a capacitance 
and a charge relaxation resistance. Their measurements 
have confirmed the prediction [Sl[7j of a quantized charge 
relaxation resistance Rq = h/2e'^ where e is the elec- 
tron charge and h the Planck constant. This prediction 
was recently shown [51 |S] to be valid for all interaction 
strength. Ref. [S] also predicted the emergence of an ad- 
ditional universal resistance Rq = h/e^ in the case of a 
large dot. 

The quantum RC circuit is described by the Ander- 
son model [lOj when the level spacing is sufficiently 
large and electron transport is not spin-polarized in con- 
trast to Refs. pH9|. In that case, the gate voltage con- 
trols the dot single-particle energy £j;(< 0). In addi- 
tion to being experimentally relevant, for the transport 
through short nanotubes small quantum dots or 

even molecules |12j . the Anderson model is fascinating 
because it displays features of strong correlation with the 
emergence of Kondo physics at low energy. The question 
of how these correlations affect the quantization of the 
charge relaxation resistance is a fundamental issue. 

The linear charge response of the quantum dot to a 
gate voltage oscillation defines the capacitance Cq and 



the resistance Rq via the low frequency expansion |5] 

^— = Co + iwCoitg + C^w ) (1) 

where h the number of electrons on the dot. The capaci- 
tance is thus the static response of the dot. The product 
LuCgRq describes relaxation towards the changing ground 
state that implies energy dissipation [51 [T3| . In this Let- 
ter, we investigate the dynamical charge response of the 
Anderson model at zero temperature and finite magnetic 
fields and we evidence a giant charge relaxation resis- 
tance phenomenon associated with the destruction of the 
Kondo effect at intermediate fields. 

More precisely, by applying a low-energy effective 
Fermi liquid theory |15[ 116) . we derive a generalized 
Korringa-Shiba formula [T7| for the charge susceptibil- 
ity that extends to finite magnetic fields. The charge 
relaxation resistance then depends only on static suscep- 
tibilities that are computed analytically resorting to the 
Bethe ansatz solution [THl dHl HI] in the Kondo regime. 
At zero magnetic field, the original Korringa-Shiba [17] 
formula predicts the quantized value Rq — h/Ae'^ inde- 
pendent of interactions. This result agrees with the non- 
interacting scattering approach with two (spin) conduct- 
ing channels [6,^7^. At large magnetic fields, the dot be- 
comes spin-polarized, reducing electron transfer in both 
spin channels, and the quantized value Rq = h/Ae^ is re- 
covered perturbatively. In the crossover regime between 
these two limits, a peak was observed in the Numerical 
Renormalization Group (NRG) calculations of Ref. [13] , 
where it is attributed to spin fluctuations in the dot. 
Hereafter, we derive analytically the emergence of this 
peak in the Kondo regime and derive its scaling proper- 
ties. In particular, the peak is found to disappear com- 
pletely at the particle- hole symmetric point. 

The origin of the peak in the resistance is related to 
the destruction of the Kondo singlet by the magnetic 
field which gives more fiexibility to the spin configura- 
tion, while the charge remains frozen by interactions. As 
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a result, a change in the gate voltage significantly mod- 
ifies the magnetization, with an increase of dissipation 
by particle-hole excitations. The charge however is rela- 
tively insensitive to the gate voltage and the capacitance 
remains small. An increasing dissipation oc G^Rq with 
an almost constant capacitance Co thus leads to an in- 
creasing charge relaxation resistance Rq. A further in- 
crease of the magnetic field eventually polarizes the spin 
on the dot, reduces spin flexibility and thereby energy 
dissipation. Hence, the charge relaxation resistance Rq 
passes through a maximum when the Zeeman energy is 
comparable to the Kondo energy. 

The Hamiltonian of the Anderson model is given by 

cr,k <y , ^ 

with ho- = d\,dcr the number of spin a electron on the dot, 
fi ~ fi^ + fi^, the linear spectrum Ska = — go'^J-BB/2 
of conduction electrons characterized by the constant 
density of states vq, and the energy levels Eda — £d ^ 
ga^BB/2 of the dot. /is is the Bohr magneton, g the 
Lande factor, B the external magnetic field and a — ± 
refers to fj i states, respectively. The two terms in the 
second line of Eq. ^ describe respectively Coulomb in- 
teraction and tunneling from the dot to the lead with the 
hybridization constant F — Trvt^t^ . In the presence of an 
AC drive of very small amplitude, £d = e^+Scj cos with 
Euj ^ the system relaxes towards the evolving ground 
state of the Hamiltonian and the dissipated power 

V =]^elujlTnXc{w). (3) 

which is given by linear response theory, is proportional 
to the imaginary part of the dynamical charge suscepti- 
bility xc{t - t') = ie{t - t'){[h{t), h{t')]). 

NRG calculations [16] and RG arguments [20] have 
shown that the low energy properties of the Anderson 
model ([2]) are always those of a Fermi liquid. The effec- 
tive Fermi liquid Hamiltonian takes the form 

H = ^e-koa\oa'ka + ^ Ka{'£d)a\o°'k'<y (4) 

cr,k k,k' ,cr 

The free quasiparticles of the first term are related to 
the original fermions Cka by a phase shift of 7r/2. The 
second term is a marginal perturbation corresponding to 
a potential scattering at the impurity site. It defines a 
line of fixed points parametrized by £d connecting the 
Kondo regime (for £d — —U/2) to the mixed valence 
regimes (for — or ~ — J7). The potential is related 
to the mean occupation of the dot via the Friedel sum 
rule {ho) =1/2— ^ aictm\['KVQKo{£d)\- Note that the 
potentials Ko{£d) formally also depend on U , F and B. 



Again we study the response to the AC drive = + 
EcjCoswt with EujtOJ 0. Expanding the potentials as 
Ka{£d) = Ko + K'oie^d) cosujt, we change the basis to 
the one-particle states [TB] that diagonalize the potential 
scattering terms = Ko{e^). The remaining scattering 
term in the Hamiltonian is given by 

<7 ^ ' k,k' 

with the new quasiparticles a,k'a- The derivative of the 
occupation numbers with respect to Ed in the Friedel sum 
rule formula above introduces the static spin-dependent 
susceptibilities Xca = —d{ho)/dEd- Once inserted into 
Eq. ([5|, we obtain 

— — — 

cr k.k' 

In the static case = 0, the second term in Eq. ^ adds 
the phase shift So — —T^vaEt^Xca/vf). The Friedel sum 
rule translates it into a shift in the occupations 5 {ho) = 
~Xca £u> in agreement with the definition of the charge 
susceptibilities. The Hamiltonian in Eq. ^ is extremely 
general and it only assumes a low-energy Fermi liquid 
fixed point. A similar model can be found in Ref. [^ 
where the spin susceptibility is discussed. 

Interestingly, the low energy model Eq. ([6| provides 
an alternative to compute the dissipated power Eq. 
Following standard linear response theory, it involves the 
operators Ao = (Xco-/t'o) Z]fc,fc' ^Lofcv, coupled to the 
AC excitation in Eq. ([6]), namely 

(7 

with XAaS^ - = *^(^ - t'){[Ao{t),Ao{t')]). The op- 
erators Ao create particle-hole pairs that are responsible 
for energy dissipation. The calculation is straightforward 
and gives, at zero temperature, Imx^^(cj) = vrXco-'^: 
proportional to the density of available particle-hole pairs 
with energy lo. An identification of Eqs. (|3| and ([?]) fi- 
nally results in our generalized Korringa-Shiba formula 

Imxc(w) = TTCJ (xcf + xli) , (8) 

obtained to lowest order |22j in uj. The physical mean- 
ing of this expression is explicit. In the presence of the 
AC driving applied to the gate voltage, relaxation is nec- 
essary to adjust the occupation numbers to the instant 
ground state of the Hamiltonian. This relaxation is real- 
ized by particle-hole excitations, in each spin sector inde- 
pendently, with amplitudes (see Eq. ([6|) that are deter- 
mined by the static charge susceptibilities Xca controlling 
the variations of the spin populations with the gate volt- 
age. Eq. ([8| simply states that the energy dissipated 
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in the relaxation mechanism increases quadratically with 
these amplitudes as a result of the Fermi golden rule. 

At a general level, the low frequency properties of 
the quantum RC circuit characterized by Eq. ([T]) derive 
from the knowledge of the dynamical charge suscepti- 
bility since {h{uj)) = — Xc('^)ed(w) in the linear regime. 
The capacitance Cq = e^Xc is solely determined by the 
static charge susceptibility Xc = Xc(w -J> 0) = Xct + Xci 
that is calculated using Bethe ansatz. Hence, measuring 
the capacitance realizes a charge spectroscopy |23j. At 
zero magnetic field and large enough interaction, [/ > F, 
Xc develops a double-peak structure as a function of Ed- 
maximum in the valence regimes with a valley in the 
intermediate Kondo regime around Sd = —U/2 where 
Xc — Sr/nU^ for ?7 3> F. This strong reduction of ca- 
pacitance, or charge sensitivity, characterizes the Kondo 
limit where the charge is frozen. It contradicts the non- 
interacting scattering theory [H |6] where the capacitance 
is proportional to the density of state and would there- 
fore reveal the Kondo resonance [13] . The two approaches 
are nonetheless reconciled by noting that, the Kondo res- 
onance is mostly exhausted by spin fluctuations, and the 
density of states of charge excitations of the Anderson 
model, the holons, reproduces [21] the exact value of the 
charge susceptibility and thus of the capacitance. 

The Korringa-Shiba Eq. ([8| substituted in the expan- 
sion Eq. ([T]) expresses the resistance Rq in terms of static 
susceptibilities, computable by Bethe ansatz. At zero 
magnetic field, Xc-f — Xci = Xc/2, Eq. Q reproduces the 
standard Korringa-Shiba relation and the charge relax- 
ation resistance is found to be quantized and universal. 
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in agreement with the scattering approach involving two 
equivalent spin channels [H [?]• In the general case, we 
introduce the charge magneto-susceptibility Xm = Xct ~ 
Xci which measures the sensitivity of the magnetization 
to a change in the gate voltage. The resistance reads 



(10) 



For Ed = —U/2, particle-hole symmetry implies that the 
magnetization is extremal with respect to the gate volt- 
age and Xm identically vanishes. Eq. (|9| is thus obtained 
for all ratios of U/T. 

In the rest of this Letter, we focus on the Kondo 
regime t7 ^ F where the gate voltage explores the val- 
ley between the Coulomb peaks located around Sd — 
and £d — —U. Far enough from these Coulomb peaks, 
kdl/r ^ ln(L//F), the charge on the dot remains of order 
one and the renormalization |20j of the peak positions is 
negligible. The charge susceptibility is computed pertur- 
batively at zero magnetic field, 
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FIG. 1. (a) Scaling function $(a:) (full line) and (b) Envelope 
function F{y) (see main text). F{0) = at the particle- 
hole symmetric point y = where dTK/dsd = 0. y = ±1 
correspond to the two Coulomb peaks in transport or valence 
regime. Green circles {ed = —0.15) and blue crosses (sd = 
-0.1) are extracted from Fig. 3a) of Ref. |T3j, with U = 0.4 
and r = 0.02, by implementing Eq. (14 1 and rescaling the 
X axis, (c) Charge relaxation resistance for SdjU = —1/2 ± 
0.1967 and various ratios of U jV = 15, 10, 7. 



and remains constant as the magnetic field is increased 

with gp^sB < Vkdir. 

For moderate magnetic fields gusB <C i/fedlT, the 
magnetization of the dot is known [TOKIH] from the Bethe 
ansatz solution of the Anderson model. In the Kondo 
limit, it exhibits the scaling form, 



\ kbTk J 



(12) 



where Tk = 2^UT /iree-x^l-Kedied + C/)/2J7F] is the 
Kondo temperature and the scaling function f{x) con- 
nects the asymptotes f{x) — x/^/2'Ke for x ^ 1, and 
f{x) = 1/2 — l/(41nx) for a; 1, i.e., for low and large 
magnetic fields, e referring to Euler's number. The de- 



pendence of the magnetization Eq. ( 12 ) on the gate volt- 
age, or Ed, is via the Kondo temperature. Computing the 
derivative of the Kondo temperature with respect to Ed, 
then one finds 
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(13) 



The charge magneto-susceptibility Xm is a an odd func- 
tion of Ed + U /2 that vanishes at the particle-hole sym- 
metric point. The scaling function ^{x) = xf'{x) is rep- 
resented Fig. [ija) in good agreement with Ref. [13] . It 
exhibits a peak at = 1.0697 with $(a;o) = 0.1257. 
Inserting the results Eqs. (11) and (13 1 into Eq. (10 1, 



the scaling form of the charge relaxation resistance is ob- 
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tained, 



U 



F{y) mx)Y 



(14) 



with y = [2ed + U)/U, x = g^MsB /ksTK- The peak in 
the resistance as a function of the magnetic field is de- 
scribed by the scahng function $(a;). It is also weighted 
by the envelope F{y) = (7rV8)2y2(y2 _ i^Ajf^-^ _^ ^^2^2^ 

shown Fig. [ijb). The agreement with Ref. [13], where 
U /T — 20 is finite, is here only approximate. The 
global maximum in the resistance is thus obtained for 
ea/U = -1/2 ± 0.1967, gtisB ^ 1.0697 /csTi^, with 
Rq = 0.00142 {h/4e'^){U/ry which predicts a strong in- 
crease of the resistance maximum with the ratio U/T, as 
seen Fig. [ijc). 

For large magnetic fields g^isB '3> y^^\ed\r, the free 
orbital regime [50] is reached and straightforward per- 
turbation theory applies. The result is 
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1 r 
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with Em = This leads to Xm = for = —U/2 

as expected and, for very large magnetic fields gusB 3> 
(|£ci|,ed + U), the quantized value Eq. ^ is recovered 
for all gate voltages. Note that the standard result [BHS] 
Rq = h/2e^ is only recovered for a fully polarized Fermi 
sea in the lead. 

To summarize, the peak in the charge relaxation resis- 
tance is due to the enhancement of Xm while the total 
charge remains quenched and Xc small. In the presence 
of a finite magnetic field, the Kondo state is a mixture 
of singlet and triplet spin configurations controlled by 
the ratio of the Zeeman energy to the Kondo energy. A 
change in the gate voltage modifies this ratio and, while 
keeping the total charge of the dot almost constant, redis- 
tributes the spin up and down occupations. This leads 
to a larger number of particle-hole excitations for each 
spin species and therefore increases dissipation. At the 
particle-hole symmetric point, the Kondo energy is sta- 
tionary with respect to the gate voltage such that no spin 
redistribution occurs and the peak in the resistance is ab- 
sent. It is worth mentioning that the predicted peak in 
the charge relaxation resistance occurring at intermediate 
magnetic fields can be observed using current technology 
[3]. This work can be extended in various directions, by 
considering either Zeeman effects on a large cavity char- 
acterized by several energy levels [25] or a large number 



of channels [26]. We finally stress that our result ( 10 ) for [23] 



the resistance is valid not only in the Kondo regime but 
for all values oi U, ej, and B. 
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